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Abstract Prominences or filaments are cool clouds of partially ionized plasma living 
in the solar corona. Ground- and space-based observations have confirmed the pres- 
ence of oscillatory motions in prominences and they have been interpreted in terms 
of magnetohydrodynamic (MHD) waves. Existing observational evidence points out 
that these oscillatory motions are damped in short spatial and temporal scales by 
some still not well known physical mechanism(s). Since prominences are partially 
ionized plasmas, a potential mechanism able to damp these oscillations could be 
ion-neutral collisions. Here, we will review the work done on the effects of partial 
ionization on MHD waves in prominence plasmas. 



1 Introduction 



Quiescent solar filaments are clouds of cool and dense plasma suspended against 
gravity by forces thought to be of magnetic origin. High-resolution H a observations 
(ifll, ||2l ) have revealed that the fine structure of filaments is apparently composed by 
many horizontal and thin dark threads (see J3), for a review). The measured average 
width of resolved thin threads is about 0.3 arc. sec (~ 210 km) while their length is 
between 5 and 40 arc. sec (~ 3500 - 28000 km). The fine threads of solar filaments 
seem to be partially filled with cold plasma [ 1 1, typically two orders of magnitude 
denser and cooler than the surrounding corona, and it is generally assumed that they 
outline their magnetic flux tubes ||4][5][Tl|6l|7][8]. This idea is strongly supported by 
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observations which suggest that they are inclined with respect to the filament long 
axis in a similar way to what has been found for the magnetic field ( ll9l fTOl ITTTl > . 

Small amplitude oscillations in prominences and filaments are a commonly ob- 
served phenomenon. The detected peak velocity ranges from the noise level (down 
to 0.1 km s _1 in some cases) to 2-3 km s _1 . The observed periodic signals are 
mainly detected from Doppler velocity measurements and can therefore be associ- 
ated to the transverse displacement of the fine structures lfl2l . Two-dimensional ob- 
servations of filaments |[T3l[T4i revealed that individual threads or groups of threads 
may oscillate independently with their own periods, which range between 3 and 
20 minutes. Furthermore, [15| have shown evidence about traveling waves along a 
number of filament threads with an average phase velocity of 12 km s , a wave- 
length of 4" (~ 2800 km), and oscillatory periods of the individual threads that vary 
from 3 to 9 minutes. 

Observational evidence for the damping of small amplitude oscillations in promi- 
nences can be found in [16|. Observational studies have allowed to obtain some 
characteristic spatial and time scales. Reliable values for the damping time have 
been derived, from different Doppler velocity time series by ifTTl . in prominences, 
and by [5] in filaments. The values thus obtained are usually between 1 and 4 times 
the corresponding period, and large regions of prominences/filaments display simi- 
lar damping times. 

Finally, small amplitude oscillations in quiescent filaments have been interpreted 
in terms of magnetohydrodynamic (MHD) waves lfl"8l and, in many cases, theoret- 
ical works studying the damping of prominence oscillations have studied first the 
effect of a given damping mechanism on MHD waves in a simple, uniform, and 
unbounded media before to introduce structuring and non-uniformity. This is the 
approach that we will follow in this paper. 



2 MHD waves in unbounded partially ionized prominence 
plasmas 

Since the temperature of prominences is typically of the order of 10 4 K, the promi- 
nence plasma is only partially ionized. The exact ionization degree of prominences 
is unknown and the reported ratio of electron density to neutral hydrogen density 
lfl9ll covers about two orders of magnitude (0.1 - 10). Partial ionization brings the 
presence of neutrals in addition to electrons and ions, thus collisions between the 
different species are possible. Because of the occurrence of collisions between elec- 
trons with neutral atoms and ions, and more importantly between ions and neutrals, 
Joule dissipation is enhanced when compared with the fully ionized case. A partially 
ionized plasma can be represented as a single-fluid in the strong coupling approxi- 
mation, which is valid when the ion density in the plasma is low and the collision 
time between neutrals and ions is short compared with other time-scales of the prob- 
lem. Using this approximation it is possible to describe the very low frequency and 
large-scale fluid-like behaviour of plasmas l20l . 



MHD waves in partially ionized prominence plasmas 



3 



Partial ionization affects the induction equation, which contains additional terms 
due to the presence of neutrals and a non-zero resistivity ll2D . These additional 
terms account for the processes of ohmic diffusion, ambipolar diffusion, and Hall's 
magnetic diffusion. Ohmic diffusion is mainly due to electron-ion collisions and 
produces magnetic diffusion parallel to the magnetic field lines; ambipolar diffu- 
sion is mostly caused by ion-neutral collisions and Hall's effect is enhanced by 
ion-neutral collisions since they tend to decouple ions from the magnetic field while 
electrons remain able to drift with the magnetic field [22] . Due to the presence of 
neutrals, perpendicular magnetic diffusion is much more efficient than longitudinal 
magnetic diffusion in a partially ionized plasma. It is important to note that this is 
so even for a small relative density of neutrals. 



2.1 Homogeneous and unbounded prominence medium 

Several studies have considered the damping of MHD waves in partially ionized 
plasmas of the solar atmosphere |23l l24l |25l [26). In the context of solar promi- 
nences, lETl derived the full set of MHD equations for a partially ionized, one-fluid 
hydrogen plasma and applied them to the study of the time damping of linear, adi- 
abatic fast and slow magnetoacoustic waves in an unbounded prominence medium. 
This study was later extended to the non-adiabatic case, including thermal conduc- 
tion by neutrals and electrons and radiative losses 11271 . ETI considered a uniform 
and unbounded prominence plasma and found that ion-neutral collisions are more 
important for fast waves, for which the ratio of the damping time to the period is 
in the range 1 to 10 5 , than for slow waves, for which values between 10 4 and 10 8 
are obtained. Fast waves are efficiently damped for moderate values of the ioniza- 
tion fraction, while in a nearly fully ionized plasma, the small amount of neutrals is 
insufficient to damp the perturbations. 

In the above studies, a hydrogen plasma was considered, although 90% of the 
prominence chemical composition is hydrogen while the remaining 10% is he- 
lium. The effect of including helium in the model of l27l was assessed by ll28ll . 
The species present in the medium are electrons, protons, neutral hydrogen, neu- 
tral helium (He I) and singly ionized helium (He II), while the presence of He III is 
neglected 11291 . 

The hydrogen ionization degree is characterized by which is equivalent to the 
mean atomic weight of a pure hydrogen plasma and ranges between 0.5 for fully 
ionized hydrogen and 1 for fully neutral hydrogen. The helium ionization degree is 
characterized by 8n e — y^-, where ^Heii an d ^Hei denote the relative densities of 

SHel 

single ionized and neutral helium, respectively. Figure Q] a, b, c displays Vd/P as a 
function of the wavenumber, k, for the Alfven, fast and slow waves. In this Figure, 
the results corresponding to several helium abundances are compared for hydrogen 
and helium ionization degrees of /2h = 0.8 and 5n e = 0.1, respectively, and it can 
be observed that the presence of helium has a minor effect on the results. 
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Fig. 1 Wave damping by ion-neutral effects in a uniform medium, (a)-(c) Ratio of the damping 
time to the period, Tb / P, versus the wavenumber, k, corresponding to the Alfven wave, fast wave 
and slow wave, respectively, (d) Damping time, T/j, of the thermal wave versus the wavenumber, 
k. The different linestyles represent the following abundances: £uel = 0% (solid line), ^Heii = 10% 
(dotted line) and ^Hei = 20% (dashed line). In all computations, p.^ = 0.8 and 5n e = 0.1. The 
results for ^Hei = 10% and Sue = 0.5 are plotted by means of symbols for comparison. The shaded 
regions correspond to the range of typically observed wavelengths of prominence oscillations. In 
all the figures shown, the angle, 0, between the wavevector and the magnetic field is 7t/4. From 
(28) 



The thermal mode is a purely damped, non-propagating disturbance (a), = 0), 
so only the damping time, %d, is plotted (Figure [TJl). We observe that the effect of 
helium is different in two ranges of k. For k > 10~ 4 m -1 , thermal conduction is the 
dominant damping mechanism, so the larger the amount of helium, the shorter Tq 
because of the enhanced thermal conduction by neutral helium atoms. On the other 
hand, radiative losses are more relevant for k < 10~ 4 m~ 1 . In this region, the thermal 
mode damping time grows as the helium abundance increases. Since these variations 
in the damping time are very small, we again conclude that the damping time ob- 
tained in the absence of helium does not significantly change when helium is taken 
into account. Therefore, the inclusion of neutral or single ionized helium in partially 
ionized prominence plasmas does not modify the behaviour of linear, adiabatic or 
non-adiabatic MHD waves already found by [21 1 and 11271 . On the other hand, in 
Figure QJ we can observe that in the case of slow waves, and within most of the 
interval of observed wavelengths in prominence oscillations, the ratio between the 
damping time and the period agrees with the observational determinations, which is 
due to the joint effect of ion-neutral collisions and non-adiabatic effects ([27], [28 1). 
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Fig. 2 Sketch of the longitudinally homogeneous prominence thread model used in 1 30 1 with / = 
and in 13111321 with / ^ 0. In this Figure a denotes the radius of the cylinder while in the text we 
use R. Adapted from 1 321- 



3 Magnetohydrodynamic Waves in Prominence Threads 

In this Section, we summarize the results of papers which investigate the damping 
of MHD waves in partially ionized prominence thread models . For simplicity, early 
investigations neglected the variation of density along the thread and took into ac- 
count the variation of density in the transverse direction only. Subsequent works 
incorporated longitudinal inhomogeneity in addition to transverse inhomogeneity. 



3.1 Longitudinally Homogeneous Thread Models 

The first papers that studied partial ionization effects on wave propagation in a lon- 
gitudinally homogeneous prominence thread model were [30, 31 . 32']. These authors 
investigated linear MHD waves superimposed on a straight magnetic cylinder of ra- 
dius R, representing the thread itself, and embedded in a fully ionized and uniform 
coronal plasma. Gravity was neglected and the magnetic field was taken constant 
along the axis of the cylinder. ||30l considered an abrupt jump of density in the 
transverse direction from the internal (prominence), p p , to the external (coronal), 
p c , densities at the thread boundary, while OTl l32l replaced the discontinuity in 
density by a continuous variation of density in a region of thickness /. For / = the 
equilibrium of ||3Tl[32l reverts to that of ||30l . Hence the ratio l/R indicates the in- 
homogeneity length-scale in the transverse direction. In both papers the prominence 
plasma was assumed partially ionized with an arbitrary ionization degree, while 
the external coronal medium was fully ionized. The single-fluid approximation was 
adopted and Ohm's, Hall's, and Cowling's terms were included in the induction 
equation. Thus, the equilibrium configuration is similar to the classical straight flux 
tube model investigated by, e.g., [33 34], with the addition of partial ionization. A 
sketch of the model is displayed in Figure[2] 
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In this model the observed transverse oscillations of prominence threads can be 
interpreted in terms of transverse (Alfvenic) kink modes. Because of their observa- 
tional relevance, here we discuss the results for transverse (Alfvenic) kink waves 
only. The interested reader is refereed to the original papers 11301 [3T1 [32l where the 
results of other waves are explained in detail. It is well known that for / ^ the kink 
mode is resonantly coupled to Alfven continuum modes in the region of transversely 
non-uniform density. As a consequence the kink mode is damped by resonant ab- 
sorption. In addition, the kink mode is also damped by magnetic diffusion effects 
due to partial ionization. In the fully ionized case, the ideal resonant damping of the 
kink mode in prominence threads was investigated by l35l[36ll . while partial ioniza- 
tion does not affect the mechanism of resonant absorption, which is an ideal process 
independent of dissipation by ion-neutral collisions. This has been shown by ||37l 
using multifluid theory. 

11301 [31 1 studied temporal damping of standing waves. By neglecting the effects 
of Ohm's and Hall's diffusion in comparison to that of Cowling's diffusion, approx- 
imate expressions for the period, P, and for the ratio of the damping time, Td, to the 
period of the kink mode can be obtained in the long-wavelength limit, i.e., X/R^S> 1, 
where X is the wavelength. The long-wavelength limit is a reasonable approximation 
since wavelengths typically observed in prominences are roughly between 10 3 km 
and 10 5 km (see 1 18 1) while the observed widths of the threads are between 100 km 
and 600 km (see [ 3 1). The expressions for P and Td/P are 



(1) 
(2) 

where va is the prominence Alfven velocity, £ = p p / p c is the density contrast, «k = 
2k jP is the kink mode frequency (with P given by Equation ([U), | n is the fraction 
of neutrals, and Vi n the ion-neutral collision frequency. ^ n = for a fully ionized 
plasma and £, n = 1 for a neutral medium. To perform a check, we take X = 10 4 km, 
va = 50 km , and £ = 200. Equation (Q~|i gives P « 2.4 min, which is consistent 
with the observed periods. 

Regarding damping, the first term within the parenthesis of Equation (f2]i is due 
to resonant absorption and the second term is due to Cowling's diffusion. Note that 
the original expression of Tq/P given in [ 3 1 1 involves Cowling's diffusivity, 77c- In 
the present discussion we have replaced rjc by its expression in terms of <^ n and 
Vi n ( see the expression of 77c in, e.g., 1 30 ] ). Our purpose is to show that the term 
related to Cowling's diffusion is proportional to the ratio ©k/ V; n . To perform a sim- 
ple order-of-magnitude estimation of the importance of Cowling's diffusion, let us 
take a period of 3 min and compute Vj n using <^ n = 0.5, a prominence density of 
5 x 10~ n kg m~ 3 , and a prominence temperature of 8,000 K (see the expression of 
Vj n m [ 30]). The result is co^/v m « 2.38 x 10~ 4 . This estimation indicates that the 
effect of Cowling's diffusion is negligible unless the prominence is almost neutral, 
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Fig. 3 (a) Td/P vs. k z a for the kink mode in a longitudinally homogeneous thread with l/a = 
(dotted), l/a = 0.1 (dashed), l/a = 0.2 (solid), and l/a = 0.4 (dash-dotted). The shaded zone de- 
notes realistic wavelengths. Adapted from 1311 . (b) Td/P vs. Lp/L for the kink mode in a longitu- 
dinally inhomogeneous thread with l/a = 0.05 (dotted), l/a = 0.1 (dashed), l/a = 0.2 (solid), and 
l/a = 0.4 (dash-dotted). Symbols are the result from Equation {5J with l/a = 0.2. In this Figure a 
denotes the radius of the cylinder while in the text we use R. Adapted from 1 38 1. 



i.e., §n w 1, which is an unrealistic limit. Therefore, resonant absorption dominates 
the kink mode damping and the second term within the parenthesis of Equation © 
can be dropped. Hence, for £ = 200 and / /R = 0.2 we obtain Td /P ~ 3.22, which 
again is consistent with the observations. 

|[3~T1 also obtained results beyond the long-wavelength approximation by means 
of full numerical solutions. They included Ohm's and Hall's terms along with Cowl- 
ing's diffusion. iFJTI computed the kink mode Tq/P as a function of the param- 
eter k z R, where k z =2n/X (see Fig. [3^). They concluded that Hall's term is al- 
ways negligible in prominence conditions, Ohm's diffusion is only important for 
extremely long wavelengths (very small values of k z R), and Cowling's diffusion 
is only relevant for short wavelengths (large k z R). For realistic wavelengths, i.e., 
10~ 3 < k z R < 10 _1 , resonant absorption determines the damping rate of the kink 
mode and the analytical formula given in Equation (O is very accurate. 

Subsequently, ||32l used the same model as ||3T| to study spatial damping of kink 
waves. The results of [32 1 are qualitatively equivalent to those of 11311 , i.e., resonant 
damping dominates for realistic frequencies whereas Cowling's diffusion is efficient 
for high frequencies only. 



3.2 Longitudinally Inhomogeneous Thread Models 

A longitudinally homogeneous cylinder is a crude representation of a prominence 
fine structure. High-resolution observations of prominences (see [3|) suggest that 
the cool and dense material of the threads only occupies a small part of longer mag- 
netic flux tubes, with the rest of the magnetic tube probably filled with hot coronal 
plasma. The observed length of the threads is believed to be only a small percentage 
of the total length of the magnetic tube, whose feet are rooted in the solar photo- 
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Fig. 4 Sketch of the longitudinally inhomogeneous prominence thread model used in 1 38 1. In this 
Figure a denotes the radius of the cylinder while in the text we use R. Adapted from 1381 . 



sphere |[39l l40l . This observational evidence may be omitted to study propagating 
waves in the dense part of the magnetic tube if the wavelengths are much shorter 
than the length of the threads. For this case the longitudinally homogeneous mod- 
els discussed in Section [3TTI may be appropriate. However, in the case of standing 
modes, the associated wavelengths are of the order of the total length of magnetic 
field lines. Therefore the longitudinal structuring of the prominence magnetic tube 
cannot be neglected when standing modes are investigated. 

This observational evidence has been taken into account in some works which 
studied ideal, undamped kink modes in the fully ionized case (see, e.g., [41 1 421 1431 ). 
The first paper that incorporated the effects of damping by Cowling's diffusion 
and resonant absorption on standing kink modes in longitudinally inhomogeneous 
threads was l38l . These authors used the model displayed in Figure [4] It is com- 
posed of a straight magnetic cylinder of length L and radius R whose ends are fixed 
at two rigid walls representing the solar photosphere. The magnetic field is constant. 
The cylinder is composed of a region of length L p and density p p , representing the 
prominence thread, surrounded by two regions of density p e representing the evac- 
uated part of the tube. The external coronal density is p c and for simplicity it is set 
p e = p c - The prominence thread is transversely inhomogeneous in a region of thick- 
ness I where the density continuously varies from p p to p c . The prominence plasma 
is partially ionized while the coronal and evacuated plasmas are fully ionized. 

To investigate standing kink modes analytically, [38] used the thin tube approx- 
imation, i.e., R/L <C 1 and R/L p <C 1. To check this approximation we take the 
values of R and L p typically reported from the observations (see Q) and assume 
L ~ 10 s km, so that R/L v and R/L are in the ranges 2 x 10~ 3 < R/L p < 0.1 and 
5 x 10~ 4 < R/L < 3 x 10~ 3 , meaning that the use of the TT approximation is jus- 
tified. [44] have shown that the results of [38 1 remain valid beyond the thin tube 
approximation. [ 38 1 derived approximate expressions for P and the ratio Tq/P. The 
expression for P is 
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where here va is the Alfven velocity of the dense, prominence plasma only, and 
£ = p p / p c is the density constrast as before. In Equation ([3]) it is assumed that the 
prominence thread is located at the center of the magnetic tube (a general expression 
is given in [38]). A direct comparison of Equations ((TJ and © shows that the effect 
of the longitudinal structuring of the tube is to select a particular value of the wave- 
length, X, which depends of the relation between L and Lp. Regarding the damping 
rate, the expression for Td/P is not explicitly given here because it is exactly the 
same as that in Equation ©, where now (»\ l = 2k/P has to be computed using the 
period from Equation (01. Thus, as happens for kink modes in longitudinally homo- 
geneous threads, the effect of Cowling's diffusion is negligible for realistic values 
of the period when compared to that of resonant absorption. 

Figure [3p shows the ratio Td/P numerically computed by [38] as a function of 
Lp/L. Remarkably, the damping ratio is independent of the length of the thread and 
only depends on the transverse non-uniformity length scale / /R. These means that 
the expression of Tfo/P for longitudinally homogeneous (Equation (O) tubes also 
applies when longitudinal structuring is included. 
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